Spatial and Spatio-temporal Epidemiology

Urban air pollution remains a significant public health concern, contributing to premature deaths and adverse health outcomes. However, there is little causal research evaluating the effectiveness of policies designed to improve air quality. This study assesses the impact of all three stages of London's Ultra Low Emission Zone (ULEZ) on air pollution, via PM2.5 levels, and respiratory health, via prescription records for bronchodilator and respiratory corticosteroid medications. Analyses are at general practice level, using a generalised synthetic control method to estimate causal impacts. Stage 1 was associated with a statistically significant but negligible 0.77% reduction in PM2.5 levels, with no corresponding change in prescribing. Stage 2 produced a paradoxical 2.69% increase in PM2.5, alongside a 4.44% decrease in inhaled corticosteroid quantity but a 12.51% increase in average daily quantity (ADQ) usage, suggesting a worsening of disease severity among existing patients. Stage 3 yielded a 2.69% PM2.5 reduction and a modest 2.18% decrease in bronchodilator ADQ usage. Spillover effects beyond the ULEZ boundary were statistically significant, but negligible. We find overall that the ULEZ had minimal effects on both air quality and respiratory prescribing across all three stages. These findings provide new insights into the effectiveness of ULEZ policies in reducing air pollution and its associated health impacts, suggesting the zone's effects are considerably smaller than previously reported, and that integration with broader policy measures may be necessary to achieve meaningful public health gains.

Background Dengue is a major public health challenge, and predictive models are crucial for early warning systems. However, many current modeling practices rely exclusively on climatic factors or employ complex algorithms that lack the interpretability needed for informed public health decision-making. To address these shortcomings, we developed and validated a multidimensional, interpretable statistical model to predict monthly dengue incidence. Methodology/Principal Findings We used a Generalized Linear Mixed Model (GLMM) with a Negative Binomial distribution to analyze 14 years (2010-2023) of spatiotemporal data from 37 municipalities in Huila, Colombia, an endemic region. The model integrates non-linear and lagged effects of climatic, demographic, and socioeconomic factors. The final model underwent rigorous external validation on an independent test set (2021-2023). Our model demonstrated high predictive discrimination (R2 = 0.743, Spearman's {rho} = 0.657), accurately capturing the timing of epidemic outbreaks. Key findings include the identification of an optimal thermal window for transmission at 27-28{degrees}C, a threshold effect for precipitation above 800 mm, and a saturation dynamic in outbreak autocorrelation. Conclusions/Significance This mechanistically-informed statistical approach provides a robust and transparent tool for epidemiological surveillance, successfully balancing high predictive performance with the explanatory power needed for effective, data-driven public health interventions.

Motivation: With the scaling frequency and intensity of extreme heat events across the globe, it is critical for public institutions to develop early detection systems and continuous monitoring of these events and their impacts. In Brazil, Rio de Janeiro was the first city to publish its heat protocol, with the Rio Heat Dashboard as a central component of this system. Implementation: The dashboard was implemented using R/Shiny and integrates climatic and health data from multiple sources. General features: The application comprises real-time heat exposure monitoring and automatic alert level classification, which is monitored daily by multiple municipal actors and supports activation of actions specified in the heat protocol. It also features a health impact module, which lists each heat event and its impact on mortality, and primary care and emergency visits. Availability: The source for full reproducibility is available through https://github.com/joaohmorais/RioHeatDashboard.

Communities exposed to flooding can experience markedly different mental health outcomes, yet conventional resilience indicators capture only part of the social and contextual conditions that may explain this variation. This study develops a multilevel and predictive framework for examining community resilience and depressive symptoms following flood exposure in Indonesia. Data were drawn from 20,303 respondents aged 15 years and older nested within 312 communities in the Indonesia Family Life Survey (IFLS-5). Depressive symptoms were assessed using the 10-item Centre for Epidemiologic Studies Depression Scale (CES-D-10), with Rasch Partial Credit Model calibration used to examine measurement properties. Bayesian multilevel models quantified between-community heterogeneity and assessed how far observable structural resources accounted for this variation. Community resilience was represented through two complementary constructs: structural resilience, based on observable socioeconomic and social-capital resources, and Latent Community Protective Capacity (LCPC), a model-derived proxy for residual contextual variation in depressive-symptom risk. Approximately 6 percent of variation was attributable to between-community differences, while observable structural resources explained only part of this heterogeneity. Structural resilience and LCPC were weakly correlated (r = 0.155). Moderation analyses provided no clear evidence that structural resilience altered the flood-depression association, while LCPC showed a directionally consistent but uncertain buffering pattern. Predictive models incorporating community-level information improved discrimination, with the best-performing model reaching an ROC-AUC of approximately 0.71. The findings suggest that observable resource-based indices provide an incomplete account of community-level mental health vulnerability and that residual contextual measures may provide complementary information, while requiring cautious interpretation and independent validation.

Background Childhood respiratory mortality in China has fallen by over 90% in three decades alongside sustained national warming, yet national long-run evidence on temperature and child respiratory mortality is lacking. Methods We linked Global Burden of Disease (GBD) 2021 mortality estimates for China - lower respiratory infections (LRI), ages 0-19, and asthma, ages 0-24, 1990-2021 - with C-LSAT 0.5 deg gridded temperature data (1990-2019), aggregated nationally and to five climate zones. Four annual indicators (mean temperature, diurnal temperature range, seasonal amplitude, interannual variability) entered regressions of log mortality rates with Newey-West standard errors. A bootstrapped (500 resamples) quadratic model probed the minimum mortality temperature (MMT), with PM2.5-adjusted analyses and future-exposure, permutation, and detrended falsification tests. Results LRI deaths fell by 96.3% (330,194 in 1990 to 12,098 in 2021; 95% uncertainty interval 9,669-14,891) and asthma deaths by 94.9% (3,287 to 167), while mean temperature rose 0.364 deg C per decade and diurnal temperature range narrowed 0.092 deg C per decade. Baseline coefficients were large (mean temperature -1.696, SE 0.174; diurnal temperature range +2.408, SE 0.336; seasonal amplitude -0.162, SE 0.082; interannual variability +2.924, SE 1.514, per 1 deg C in log rate), but the future-exposure test failed and detrending nullified every coefficient: the associations are trend-level, and short-cycle causal effects are not identifiable. Nor was the national MMT identifiable - observed temperature support spans only 6.66-8.13 deg C, and the nominal turning point of 35.84 deg C is an extrapolation artifact (quadratic term p = 0.963). Within the observed range, warming and declining mortality moved in the same direction. Conclusions The 96% decline in childhood respiratory mortality cannot be attributed to warming. China sits on the low-temperature side of the optimum, and the marginal direction of future warming requires stronger designs to establish. The falsification framework offers a discipline for climate-health inference in China.

Inappropriate antibiotic use presents a major global health challenge, particularly in low-resource settings where access to quality care is limited but antibiotics remain relatively unrestricted. This study estimates the causal effect of frontline primary care quality on inappropriate community antibiotic use, combining detailed community-based data from approximately 100 rural villages in rural China with an instrumental variable (IV) approach embedded within a double/debiased machine learning (DML) framework. We linked objective measures of village doctor clinical practice quality, measured through unannounced standardized patient visits, to household-level antibiotic use data collected from the same villages. To identify the causal effect, we constructed multiple candidate instruments from extensive provider characteristics and used an ensemble of machine learning algorithms within a flexible DML-IV framework to approximate an optimal instrument, addressing a many-weak-instruments problem. We found that improving village provider clinical practice quality reduced both antibiotic receipt during healthcare encounters for common diseases and household antibiotic storage for future self-medication. Our findings suggest that strengthening frontline primary care quality can meaningfully reduce inappropriate community antibiotic use without restricting access to essential treatment. More broadly, this study illustrates how causal machine learning can strengthen conventional causal estimation in complex observational settings in global health economics research.

Mosquito-borne diseases pose a growing public health challenge as climate change reshapes vector population dynamics. West Nile virus (WNV), transmitted between birds and Culex mosquitoes, disproportionately affects Maricopa County, Arizona, one of the nation's highest-burden counties, yet whether models that include weather and avian dynamics improve forecast accuracy remains unclear. Using a 15-year weekly time series of mosquito abundance, mosquito infection prevalence, and human cases, we developed four mechanistic model configurations of varying complexity, from mosquito-human dynamics alone to full models incorporating avian dynamics and weather forcing. We fitted each model to the weekly-observed data, generated probabilistic 1- and 2-week-ahead forecast horizons, and evaluated forecasts against a historical baseline. All configurations fit the data equally regardless of weather or avian dynamics. However, models incorporating both birds and weather created more accurate forecasts of mosquito abundance and mosquito infection prevalence, and all configurations outperformed the baseline for forecasting human cases. Forecast accuracy was highest in summer and fall, and ensemble aggregation sometimes outperformed every individual model, stabilizing predictions across the 15-year record. These findings indicate that avian and weather dynamics are most critical for predicting mosquito-specific data, positioning this framework as a scalable tool for public health planning for WNV surveillance under climate change.

Wastewater-based epidemiology (WBE) has emerged as a powerful population-level surveillance tool, but its coverage is structurally concentrated in in-network urban areas, potentially leaving rural populations underrepresented. Routine human movement between sewered (in-network) and unsewered (off-network) areas may, however, cause wastewater treatment plant (WWTP) measurements to reflect infectious disease dynamics beyond sewer boundaries. We evaluated this hypothesis using daily clinical COVID-19 testing data (January 2021-April 2022) across New York State excluding New York City (NYC). We disaggregated weekly cases and tests into in-network (WWTP catchment area) and off-network (outside WWTP catchment area) components applied to two geographic frameworks: administrative counties (N = 53 mixed-coverage) and mobility-defined communities identified through Walktrap community detection applied to census tract-level movement networks (N = 32 mixed-coverage). In/off-network COVID-19 trends were strongly correlated under both frameworks. County-level statewide aggregate correlations were high (incidence r = 0.994, positivity r = 0.996), as were individual county correlations (median r = 0.909 and 0.932, respectively). Mobility-defined community-level statewide correlations were similarly strong (r = 0.990 and 0.992), with comparable unit-level medians (r = 0.877 and 0.894). The mobility-defined community framework provided better population balance between in-network and off-network strata (87.5% vs. 69.8% in balanced range) and a higher floor on representativeness (minimum r = 0.440 vs. 0.177). Population size was the dominant predictor of in-network/off-network alignment at both scales; wastewater infrastructure density and off-network signal variability provided additional explanatory power at the mobility-defined community level. WWTPs broadly represent COVID-19 dynamics in surrounding off-network populations, supporting their use as sentinel surveillance sites. Representativeness weakens in smaller, more rural communities, and mobility-defined communities provide a complementary framework for identifying where this occurs.

Objectives To analyse the long-term effects of a cruciate ligament (CL) injury on health and socioeconomic outcomes. Methods We used a comprehensive national injury insurance database to identify CL injuries occurring in New Zealand between 2009 and 2022, and employed a doubly robust staggered difference-in-differences research design to identify the effects of these injuries on outcomes up to 10 years after injury. The outcomes of interest were healthcare use (hospitalisations, emergency department visits, medications, knee replacement surgery for osteoarthritis), associated healthcare costs, and labour market outcomes (employment rates, income, and government benefit payments). Results We identified 61 344 CL injuries for inclusion in the analysis. Over 10-year follow-up, a CL injury resulted in increased healthcare use (0.6 more hospitalizations [95%CI 0.4 to 0.7], 1.7 more days spent in hospital [95%CI 1.3 to 2.1], 0.4 more emergency department visits [95%CI 0.3 to 0.6], 2.5 more outpatient visits [95%CI 1.8 to 3.2], and 4.7 more medications dispensed [95%CI -1.8 to 11.2]) and public healthcare costs ($7 537; 95%CI 5 888 to 9 186), reduced income (-$6 060; 95%CI -11 644 to -475), and increased benefit payments ($1 152; 95%CI 542 to 1 761). Conclusion CL injuries have long-term impacts on healthcare use and socioeconomic outcomes. Strategies to reduce the incidence of CL injuries have the potential to realise large health and economic benefits.

BACKGROUND: We used anthropometric data from electronic health records (EHRs) of Swiss childrens hospitals to evaluate growth references and estimate centile curves. METHODS: We received EHRs extracted from seven Swiss childrens hospitals and analysed two samples: all children with a height, weight, body mass index (BMI), or head circumference recording, and a subsample restricted to children without diseases potentially affecting growth, weighted to represent the general population. We calculated mean z-scores based on the World Health Organization growth references adopted for Switzerland in 2011 (CH-WHO 2011) and current Swiss growth references (Swiss 2026). We estimated sex-specific centile curves in the subsample using generalised additive models for location, scale, and shape. RESULTS: We included 213,868 children with height, 448,002 with weight, 209,244 with BMI, and 67,397 with head circumference recordings. Mean z-scores in the all children sample were (CH-WHO 2011; Swiss 2026): height (0.10; -0.19), weight (0.16; -0.09), BMI (0.04; -0.07), head circumference (-0.28, -0.28); and in the subsample: height (0.34; 0.00), weight (0.27; 0.01), BMI (0.18; 0.05), and head circumference (0.04; 0.01). The 50th height, weight, BMI, and head circumference centiles of girls and boys in the subsample closely followed those of Swiss 2026, with slightly wider 3rd and 97th centiles in infancy and adolescence. CONCLUSION: Height, weight, BMI, and head circumference centiles aligned well with the Swiss 2026 growth references in Switzerland, demonstrating that hospital EHRs could contribute to future growth references.

Clinical brain MRIs from pediatric health systems represent a viable resource for modeling early neurodevelopmental trajectories and studying neurodevelopmental risk in real-world populations. However, a limitation to date has been the performance of existing segmentation tools for measuring various brain phenotypes in clinical scans. In particular, many tools underperform in infant scans due to morphological and physical changes such as rapid myelination. Here, we introduce ClinSeg: a robust segmentation approach tailored to early-life clinical MRIs with variable orientation, resolution, and contrast. We leverage existing registration and synthetic data generation tools to construct a training corpus for a 3d U-Net spanning anatomical and contrast diversity, including scans with morphological abnormalities from a pediatric hospital. Validated against manual segmentations, ClinSeg outperforms existing models in infancy while matching them in childhood and adolescence. Finally, ClinSeg enables the construction of reference brain growth trajectories in 11,699 individuals from 0-21 years of age, leading to the detection of more nuanced age-related findings in clinical groups.

Background Rural-urban disparities in SARS-CoV-2 infection epidemiology remain poorly quantified and understood in Zambia despite differences in healthcare access, services and preventive interventions. This study examined the geographical distribution and associated factors of SARS-CoV-2 cases across selected rural and urban districts of Zambia. Methods A convergent mixed-methods study comprised of quantitative survey and qualitative interviews was conducted in; Ndola (Urban), Kafue (Peri-urban) and Lufwanyama (Rural). The proximate determinant framework guided variable selection and interpretation. Quantitative combined (Hospital-surveillance data with community survey), while qualitative included In-depth interviews. Participants were sampled using multistage sampling technique. Quantitative data were analysed using STATA version 17, while qualitative data were analysed thematically. Findings were integrated through triangulation. Results A total of 528 participants were included, with a median age 31 years (15-71). Overall SARS-CoV-2 positivity was 12.6%, varying across rural (16.5%), peri-urban (14.9%), and urban (9.9%) settings, though residence was not associated with infection (P<0.132). Participants aged [&ge;]49 years had significantly higher odds of infection (aOR=8.78; 95% CI:1.15-66.99), whereas secondary education (aOR=0.37; 95% CI:0.16-0.86) and hospital-based testing (aOR=0.37; 95% CI:0.15-0.92) were associated with lower odds of infection. Vaccine uptake was highest in urban areas but was not independently associated with infection. Qualitative findings revealed marked rural-urban differences in perceived susceptibility, testing access, vaccine decision-making, and adherence to preventive measures, explaining several quantitative observations. Conclusion SARS-CoV-2 infection across rural and urban settings in Zambia was influenced by demographic, behavioral, and health-system factors rather than geographic residence alone. These findings highlight the need for context-specific prevention strategies, equitable access to testing, strengthened community surveillance, and targeted risk communication to improve preparedness and response for future respiratory disease outbreaks.

Background Little evidence is available on the epidemiology of different knee injuries at a whole-of-population level. The objective of this article is to provide accurate estimates of knee injury incidence by harnessing the unique comprehensive, population-wide data of New Zealand's universal no-fault injury insurance provider, the Accident Compensation Corporation (ACC). Methods We obtained insurance claims data from ACC covering all knee injury insurance claims approved between 2015 and 2024. We calculated the number of injuries and the incidence rate per 100 000 population, by injury type, year, sex, ethnicity, and age. Results The total number of injuries increased from 184 710 (4 067 per 100 000 population) in 2015 to 244 155 (4 701 per 100 000) in 2024. The most common injuries were other/unspecified ligament sprains, contusions, and collateral ligament sprains. Ligament and cartilage injuries were more common for males than for females, while contusions were more common for females. Ligament tears and dislocations were more common in younger people (15 to 35 years of age), while cartilage injuries were more common at older ages (40 to 65 years). Discussion and Conclusions The rate of knee injuries observed in this study was higher than previously reported in other settings, probably due to broader coverage of injuries treated in primary and community care settings. A broad range of injuries were common, including those that have received less attention in the epidemiological literature to date. More research is needed on the prevention, burden, and outcomes of different knee injuries, beyond a narrow focus on cruciate ligament injuries.

Occupational ocular injuries are a preventable yet neglected public health problem, particularly in low- and middle-income countries. Maintenance workers are exposed to diverse ocular hazards daily, yet compliance with protective measures is consistently poor. A descriptive cross-sectional study was conducted among 85 maintenance workers at the Maintenance and Essential Services Organization (MESO) of Kwame Nkrumah University of Science and Technology (KNUST), Ghana, recruited through stratified convenience sampling across seven occupational sections. A structured questionnaire assessed knowledge of ocular hazards and protective equipment, attitudes toward ocular safety, and safety practices. Data were analyzed using IBM SPSS version 26 (IBM Corp., Armonk, NY, USA); chi-square and Fishers exact tests assessed associations (p < 0.05). Participants were predominantly male (84/85, 98.8%), with a mean age of 44.5 {+/-} 10.4 years. Overall knowledge was good (mean 9.40 {+/-} 1.59 out of 11), but attitude and practice scores were average (2.78 {+/-} 0.92 and 3.27 {+/-} 0.93, respectively). Most workers correctly identified goggles and face shields as protective, but only about half recognized that ordinary sunglasses and spectacles offer inadequate protection. Although 97.6% (83/85) recognized the need for ocular protection, only 7.1% (6/85) reported consistent protective eyewear use, and fewer than half (45.9%, 39/85) had received formal ocular safety training. Routine general protective equipment use was significantly associated with ocular protection use (Fishers exact test, p = 0.011). Sand and dust particles were the leading causes of injury and only 25% (5/20) of injured workers sought formal care. Workers demonstrated good knowledge but poor attitudes and practices toward ocular safety, suggesting that knowledge alone does not translate into protective behaviour even within a relatively well-resourced institutional setting. Findings suggest that limited access to task-appropriate protective eyewear may represent an important institutional barrier. Institutional PPE supply and section-specific safety training are essential to bridge this knowledge-practice gap.

Background: Inequalities exist for many health outcomes, but there is limited evidence regarding menstrual symptoms despite their importance for health and wellbeing. We aimed to investigate inequalities in menstrual symptoms according to socioeconomic position and childhood adversity. Methods: In two generations (G0 mothers and G1 offspring) from the Avon Longitudinal Study of Parents and Children (ALSPAC), a UK prospective cohort study, we examined associations of multiple indicators of socioeconomic position (SEP) and adverse childhood experiences (ACEs) with menstrual symptoms (pain, abnormal uterine bleeding, and premenstrual syndrome (PMS) measured 3-8-years post-birth in G0 and 17-21-years-old in G1), using multivariable logistic regression. Samples ranged from 4,828 to 9,335 G0 participants and 1,288 to 2,757 G1 participants depending on the exposure-outcome association. Missing data were addressed using multiple imputation and inverse probability weighting. Results: Financial difficulties were associated with greater odds of menstrual pain (G1 OR 1.41; 95% CI 1.07, 1.86: G0 OR 1.55; 95% CI 1.36, 1.76) and irregular cycles (G1 OR 1.60; 95% CI 1.12, 2.29: G0 OR 1.48; 95% CI 1.27, 1.72) in both generations, as well as with short/long cycle lengths in G0 only. Lower education and manual social class were also associated with these three menstrual symptoms in at least one generation. Conversely, higher SEP was associated with PMS in both generations. Higher cumulative ACEs were consistently associated with menstrual pain (4+ compared to none: G1 OR 2.15; 95% CI 1.48, 3.11: G0 OR 1.52; 95% CI 1.29, 1.80) and irregular cycles (G1 OR 1.92; 95% CI 1.20, 3.09: G0 OR 1.54; 95% CI 1.26, 1.87) but not cycle length. Lower parental education, financial difficulties, and cumulative ACEs were associated with heavy bleeding in G1 offspring only, whereas financial difficulties, own manual social class, and cumulative ACEs were associated with prolonged bleeding in G0 mothers only. Higher cumulative ACEs were also associated with PMS in G1 offspring only. Conclusions: We found evidence of inequalities according to socioeconomic disadvantage and childhood adversity for multiple menstrual symptoms, although some associations were only observed in one generation. Findings suggest that menstrual symptoms are disproportionately experienced by socially and socioeconomically disadvantaged women.

Background: Vitamin K prophylaxis is universally recommended for newborns to prevent life threatening vitamin K deficiency bleeding. Although not on the immunization schedule, vitamin K prophylaxis is often coadministered with hepatitis B birth dose and erythromycin ophthalmic ointment, and rising hesitancy around vaccines/preventive care may spill over into vitamin K administration. Methods: We conducted a retrospective cohort study using Truveta electronic health record data with linked mother-child dyads. Live births to mothers aged 15-49 from January 1, 2019 through June 30, 2026 were included. Vitamin K administration was defined as documentation on the birth date or following day. Logistic regression assessed sociodemographic predictors of non-receipt, and interrupted time series analysis evaluated changes after January 2026. Results: Among 1,026,375 infants, 995,628 (96.97%) had documented vitamin K administration. Non-receipt increased from an average of 2.1% during 2019-2022 to 4.3% in 2025 and 6.1% in 2026, reaching 8.10% in June 2026. Older maternal age, non-Hispanic or Latino ethnicity, Medicaid or unknown insurance, and year of delivery were associated with greater odds of non-receipt. After January 2026, there was no immediate step change, but the odds of vitamin K receipt declined an additional 10% per month (OR: 0.90; 95% CI, 0.88-0.91). Conclusions: Vitamin K non-receipt increased over the study period and accelerated after January 2026. Because vitamin K recommendations were not changed by the January vaccine schedule, this association may reflect broader impacts to confidence in newborn preventive care. Future studies should examine causal mechanisms, parental decision-making, and associated clinical outcomes.

Background: Children spend more than 90% of their time indoors, and early childhood education settings (ECEs) are an understudied, high-occupant-density indoor microenvironment where exposure to volatile organic compounds, particulate matter, and other toxicants has been documented. Limited knowledge exists on ECE-specific exposures affecting young children and how they compare to exposures in the home. Methods: This prospective, repeated-measures pilot study targeted enrollment of 44 preschool-aged children and 8 ECE staff across two geographically and sociodemographically distinct ECEs in metropolitan Atlanta, Georgia. Paired silicone wristbands, one home-designated and one ECE-designated, were exchanged between settings across three consecutive days and nights beginning at enrollment to characterize microenvironment-specific exposure. A single spot urine sample was also collected from each child. Continuous indoor air quality monitoring was conducted in two classrooms per site. Caregivers and ECE staff completed structured questionnaires assessing home and ECE environmental characteristics, child respiratory risk, and protocol feasibility and acceptability. Feasibility was evaluated using eight pre-specified indicators spanning recruitment and enrollment, wristband wear duration and loss by microenvironment, urine sample collection completeness, and survey completion by instrument and respondent group. Conclusion: This pilot will establish feasibility and acceptability parameters for a paired, multi-matrix silicone wristband protocol across home and ECE microenvironments. Findings will inform the design, sample size, and power calculations for a subsequent study testing indoor air interventions and pediatric respiratory outcomes in ECEs. Feasibility outcomes are reported in a companion manuscript.

Background Alcohol-associated injury is a leading cause of emergency department (ED) utilization in the United States and a clinically important driver of preventable morbidity across the adult lifespan. Prior surveillance research has characterized how the rate and severity of alcohol-associated injury vary by patient age, but whether the seasonal timing of injury risk is equally predictable across age groups (a question directly relevant to the timing of clinical screening intensification and public health intervention) has not been formally tested. Methods We conducted a retrospective surveillance analysis of 45,876 alcohol-associated ED visits among adults aged 18 years and older, identified from the National Electronic Injury Surveillance System (NEISS), 2019-2025 (weighted national estimate: 2,092,319 visits), using the structured Alcohol_Involved indicator introduced into NEISS case abstraction in 2019. Patients were stratified by sex and five age groups (18-24, 25-34, 35-49, 50-64, and [&ge;]65 years). Single-harmonic cosinor (Poisson) regression was used to estimate the seasonal peak day of injury risk (acrophase) for each stratum. To assess reliability, we performed leave-one-year-out jackknife resampling (seven iterations per group), case-resampling bootstrap confidence intervals (1,000 iterations), and likelihood-ratio tests of seasonal-phase interactions. Results Peak injury timing differed significantly across age groups (X^2 [8] = 2356.2, p < .0001). Adults aged 25-64 years showed a highly reproducible early-to-mid-July peak, with jackknife estimates shifting [&le;]14 days when any single study year was excluded. Adults aged [&ge;]65 years showed significant seasonal variation annually (all p < .0001, amplitude comparable to younger groups) but a pooled peak estimate that shifted by up to 100 days across jackknife iterations. Sex-stratified analyses revealed that this instability was driven entirely by females aged [&ge;]65 years (jackknife range: 332 days, peak consistently in late October through early January) rather than males aged [&ge;]65 (jackknife range: 31 days, peak consistently in early August). Hospital admission rates increased monotonically with age from 9.0% (18-24 years) to 31.8% ([&ge;]65 years). Conclusions Alcohol-associated injury follows a reproducible, calendar-stable summer seasonal pattern in adults aged 25-64 years. Among adults [&ge;]65 years, the previously reported temporal instability is concentrated in the female subgroup, whose seasonal injury risk does not converge on a fixed calendar window. These findings suggest that fixed-calendar prevention and screening strategies are well suited to working-age adults and older men, but older women may require a year-round, individually tailored approach. Keywords: Alcohol-related injury; Emergency department; Seasonality; Age factors; Sex differences; Injury surveillance; Cosinor analysis; Older adults

Wearable (accelerometer) data and machine-learning allow objective assessment of the amount of daily physical activity. However, wearable-derived human activity is subject to measurement error. No studies have corrected the dose-response association between physical activity and survival time to chronic diseases, including cardiovascular disease (CVD). The objective is to estimate the measurement error-corrected association between CVD events and multiple measures of daily duration of light and total physical activity, derived from machine-learning and conventional accelerometer-processing methods. Our method combined an accelerated failure time model, spline, and simulation-extrapolation (SIMEX). The method recovered the true dose-response non-linear association in simulated data, while the naive model failed to capture it due to substantial attenuation. Application to the UK Biobank accelerometer cohort also showed an increased protective association of total physical activity after SIMEX correction (Time Ratio [TR] = 1.56, 95% CI: 1.28-1.82 vs. TR = 1.38, 95% CI: 1.24-1.54 for SIMEX-corrected vs. uncorrected dose-response association between the 95th and 5th percentiles of total activity), with a similar increase for light physical activity. Sensitivity analysis indicates that the female population experiences a substantially larger protective association after SIMEX correction than males. Dose-response survival analysis is a widely used analytical method in physical activity epidemiology and benefits from measurement error correction.

Background In adults, cold-attributable mortality exceeds heat-attributable mortality roughly 17-fold. Child-specific evidence has begun to emerge only recently - a nationwide Brazilian case-crossover study located the minimum mortality temperature (MMT) for under-five deaths, and a 56-country survey-based analysis linked monthly temperature anomalies to under-five mortality - but no multi-country, climate-zone-resolved estimate of the childhood respiratory-infection MMT exists, and whether temperature variability is independently associated with childhood respiratory mortality at the global scale is unknown. We quantified both. Methods We combined Global Burden of Disease 2023 mortality estimates, lower respiratory infection (LRI) deaths at ages 0-19 years and asthma deaths at ages 0-24 years, 171 countries, 1990-2023 - with 0.5 deg monthly land temperature and diurnal temperature range (DTR) fields from C-LSAT/C-LDTR (1901-2023). Four exposure dimensions (annual mean, DTR, seasonal amplitude, interannual variability) entered two-way fixed-effects models with Driscoll-Kraay standard errors. A quadratic term in mean temperature located the MMT, with percentile confidence intervals from a 300-replication country-cluster bootstrap. Future-exposure leads, country-level detrending, and permutation tests assessed contemporaneous causality, applied to both the linear coefficients and the quadratic term generating the MMT; national pneumococcal conjugate vaccine (PCV3) coverage and ambient PM2.5 exposure series were added as time-varying mechanistic covariates. Results The childhood LRI MMT was 17.1 C (95% CI 14.7-19.8), the 36th percentile of the annual-temperature distribution; zone estimates were 24.7 C in tropical and 15.8 C in subtropical countries, with weak temperate and no subarctic identification. The quadratic term underpinning the MMT, however, failed both falsification checks - future temperatures reproduced the U-shape and country-level detrending erased it - so these MMT values describe a trend-level geographic pattern of the annual construct rather than a contemporaneous dose-response. Interannual temperature variability was positively associated with LRI (+0.278, 95% CI 0.102-0.454; p = 0.002) and asthma mortality (+0.836, 95% CI 0.447-1.226; p = 2.6 x 10^-5) per 1 C, but future-exposure models returned nearly identical significant coefficients and detrending erased significance, supporting only a trend-level association; adjustment for national PCV3 coverage and PM2.5 exposure left these estimates essentially unchanged. Annual mean temperature was likewise inversely associated with both outcomes at the trend level; DTR and seasonal amplitude showed no independent within-country effects. Conclusions This study provides the first multi-country, climate-zone-resolved geography of the optimal temperature for childhood respiratory survival, spanning 171 countries; because the underlying quadratic association is trend-level, the estimates are directional. The observed variability-mortality associations are trend-level signals rather than contemporaneous causal evidence; daily-scale, child-specific designs are required to determine whether short-term thermal variability affects paediatric respiratory mortality.